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Abstract. We discuss a few recent developments that are important for under- 
standing of MHD turbulence. First, MHD turbulence is not so messy as it is usually 
believed. In fact, the notion of strong non-linear coupling of compressible and incom- 
pressible motions along MHD cascade is not tenable. Alfven, slow and fast modes 
of MHD turbulence follow their own cascades and exhibit degrees of anisotropy 
consistent with theoretical expectations. Second, the fast decay of turbulence is not 
related to the compressibility of fluid. Rates of decay of compressible and incompress- 
ible motions are very similar. Third, viscosity by neutrals does not suppress MHD 
turbulence in a partially ionized gas. Instead, MHD turbulence develops magnetic 
cascade at scales below the scale at which neutrals damp ordinary hydrodynamic mo- 
tions. Forth, density statistics does not exhibit the universality that the velocity and 
magnetic field do. For instance, at small Mach numbers the density is anisotropic, 
but it gets isotropic at high Mach numbers. Fifth, the intermittency of magnetic 
field and velocity are different. Both depend on whether the measurements are done 
in local system of reference oriented along the local magnetic field or in the global 
system of reference related to the mean magnetic field. 
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1. What is MHD Turbulence? 

This short review is focused on the recently uncovered basic properties 
of MHD turbulence^. We also briefly deal with recovery of the 3D 
statistics of turbulent velocity from observations, which is a theoretical 
problem in itself. 

When do we expect fluids to be turbulent? A fluid of viscosity u 
becomes turbulent when the rate of viscous dissipation, which is ~ v/L'^ 
at the energy injection scale L, is much smaller than the energy transfer 
rate ~ Vl/L, where Vl is the velocity dispersion at the scale L. The 
ratio of the two rates is the Reynolds number Re = V^L/u. In general, 
when Re is larger than 10 — 100 the system becomes turbulent. Chaotic 
structures develop gradually as Re increases, and those with Re ~ 
10'^ are appreciably less chaotic than those with Re ~ 10*^. Observed 
features such as star forming clouds are very chaotic for Re > 10^. 
This makes it difficult to simulate realistic turbulence. The currently 

^ It is not possible to cite all the important papers in the area of MHD turbulence 
and turbulent molecular clouds. An incomplete list of the references in a recent 
review on the statistics of MHD turbulence by Cho, Lazarian & Vishniac (2003a; 
henceforth CLVOSa) includes about two hundred entries. 
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available 3D simulations containing 512 grid cells along each side can 
have Re up to ~ O(IO^) and are limited by their grid sizes. Therefore, 
it is essential to find ^^scaling laws" in order to extrapolate numerical 
calculations {Re ~ O(IO^)) to real astrophysical fluids {Re > 10^). We 
show below that even with its limited resolution, numerics is a great 
tool for testing scaling laws. 

Kolmogorov theory provides a scaling law for incompressible non- 
magnetized hydrodynamic turbulence (Kolmogorov 1941). This law 
provides a statistical relation between the relative velocity vi of fluid 
elements and their separation I, namely, vi ~ l^^^. An equivalent de- 
scription is to express spectrum E{k) as a function of wave number k 
{r^ l/l). The two descriptions are related by kE{k) ~ vf. The famous 
Kolmogorov spectrum is Eik) ~ k~^/^. The applications of Kolmogorov 
theory range from engineering research to meteorology (see Monin & 
Yaglom 1975) but its astrophysical applications are poorly justifled 
and the application of the Kolmogorov theory can lead to erroneous 
conclusions (see reviews by Lazarian et al. 2003 and Lazarian & Yan 
2003). 

Let us consider incompressible MHD turbulence first^. There have 
long been understanding that the MHD turbulence is anisotropic (e.g. 
Shebalin et al. 1983). Substantial progress has been achieved recently 
by Goldreich k. Sridhar (1995; hereafter GS95), who made an inge- 
nious prediction regarding relative motions parallel and perpendicular 
to magnetic field B for incompressible MHD turbulence. An important 
observation that leads to understanding of the GS95 scaling^ is that 
magnetic field cannot prevent mixing motions of magnetic field lines if 
the motions are perpendicular to the magnetic field. Those motions will 
cause, however, waves that will propagate along magnetic field lines. 
If that is the case, the time scale of the wave-like motions along the 
field, i.e. ~ 1\\/Va, {l\\ is the characteristic size of the perturbation along 
the magnetic field and Va = B / yJAirp is the local Alfven speed) will 
be equal to the hydrodynamic time-scale, l±/vi, where Z_l is the char- 
acteristic size of the perturbation perpendicular to the magnetic field. 
The mixing motions are hydrodynamic-like^. They obey Kolmogorov 



^ Traditionally there is insufficient interaction between researchers dealing with 
compressible and incompressible MHD turbulence. This is very unfortunate, as 
we will show later that there are many similarities between the properties of 
incompressible MHD turbulence and those of its compressible counterpart. 

^ Here we provide a more intuitive description, while a GS95 presents a more 
mathematical one. 

"* Simulations in Cho, Lazarian & Vishniac ((2002a, 2003b) that the mixing mo- 
tions are hydrodynamic up to high order. These motions according to Cho et al. 
(2003) allow efficient turbulent heat conduction. 



lazl.tex; 2/02/2008; 12:36; p. 2 



Properties of MHD Turbulence 



3 



scaling, vi (x l/'' , because incompressible turbulence is assumed. Com- 

2/3 

bining the two relations above we can get the GS95 anisotropy, Zy oc l_l 

(or A;|| oc k^/^ in terms of wave-numbers). If we interpret /|| as the 
eddy size in the direction of the local magnetic field, and lj_ as that in 
the perpendicular directions, the relation implies that smaller eddies 
are more elongated. The latter is natural as it the energy in hydrody- 
namic motions decreases with the decrease of the scale. As the result it 
gets more and more difficult for feeble hydrodynamic motions to bend 
magnetic field lines. 

GS95 predictions have been confirmed numerically (Cho & Vish- 
niac 2000; Maron & Goldreich 2001; Cho, Lazarian & Vishniac 2002a, 
hereafter CLV02a; see also CLV03a); they are in good agreement with 
observed and inferred astrophysical spectra (see CLV03a). What hap- 
pens in a compressible MHD? Does any part of GS95 model survives? 
Literature on the properties of compressible MHD is very rich (sec re- 
views by Pouquet 1999; Cho & Lazarian 2003b and references therein). 
Higdon (1984) theoretically studied density fluctuations in the inter- 
stellar MHD turbulence. Mattliacus & Brown (1988) studied nearly in- 
compressible MHD at low Mach number and Zank & Matthaeus (1993) 
extended it. In an important paper Matthaeus et al. (1996) numerically 
explored anisotropy of compressible MHD turbulence. However, those 
papers do not provide universal scalings of the GS95 type. 

The complexity of the compressible magnetized turbulence with 
magnetic field made some researchers believe that the phenomenon 
is too complex to expect any universal scalings for molecular cloud 
research. Alleged high coupling of compressible and incompressible 
motions is often quoted to justify this point of view (see discussion 
of this point below). 

In what follows we discuss the turbulence in the presence of regular 
magnetic field which is comparable to the fluctuating one. Therefore 
for most part of our discussion, we shall discuss results obtained for 
SV ~ SB/^Airp ~ So/-v/47rp, where SB is the r.m.s. strength of the 
random magnetic field. However, we would argue that our choice is 
not so restrictive as it may be seen. Indeed, at the scales where the 
velocity perturbations are much larger than the Alfven velocity, the 
dynamical importance of magnetic field is small. Therefore we expect 
that at those scales turbulent motions are close to hydrodynamic ones. 
At smaller scales where the local turbulent velocity gets smaller than 
the Alfven speed we believe that our picture will be approximately true. 
We think that the local magnetic field should act as Bq , while the small 
scale perturbations happen in respect to that local field. This reasoning 
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Fourier 
transform 




real space Fourier space Bq''^ plane 

Figure 1. Separation method. We separate Alfven, slow, and fast modes in Fourier 
space by projecting the velocity Fourier component Vk onto bases ^a, £.s, and ^f, 
respectively. Note that = —<P- Slow basis and fast basis ^/ lie in the plane 
defined by Bo and k. Slow basis lies between —9 and ky . Fast basis ^/ lies between 
k and kx. 



is in agreement with calculations in Clio, Lazarian & Visliniac (2003b) 
and Cho &; Lazarian (2003a). 



2. Does the Decay of MHD Turbulence Depend on 
Compressibility? 

Many astrophysical problems, e.g. the turbulent support of molecular 
clouds (see review by McKee 1999), critically depends on the rate of 
turbulence decay. For a long time magnetic fields were thought to be 
the means of making turbulence less dissipative. Therefore it came as 
a surprise when numerical calculations by Mac Low et al. (1998) and 
Stone, Ostriker, & Gammie (1998) indicated that compressible MHD 
turbulence decays as fast as the hydrodynamic turbulence. This gives 
rise to a erroneous belief that it is the compressibility that is responsible 
for the rapid decay of MHD turbulence. 

This point of view has been challenged in Cho & Lazarian (2002, 
2003a, henceforth CL02 and CL03, respectively). In these papers a 
technique of separating different MHD modes was developed and used 
(see Fig. 1). This allowed us to follow how the energy was redistributed 
between these modes. 

Should the different MHD modes be strongly coupled when the 
turbulence is strong? A naive answer is "yes". Indeed, strong turbu- 
lence implies strong field line wondering. This mixes up Alfven and 
fast modes. In addition, one can show through calculations that the 
magnetic non-linearities result in the drainage of energy from Alfvenic 
cascade. However, a remarkable feature of the GS95 model is that 
Alfven perturbations cascade to small scales over just one wave period, 
which gets shorter and shorter as we move along the cascade. The 
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competing effects coupling different modes usually require more time^ . 
Wc note that as the consequence of this reasoning we should assume 
that the properties of the Alfvenic cascade (incompressible cascade!) 
should not strongly depend on the sonic Mach number. 

Are the arguments above correct? The generation of compressible 
motions (i.e. radial components in Fourier space) from Alfvenic turbu- 
lence is a measure of mode coupling. How much energy in compressible 
motions is drained from Alfvenic cascade? According to closure calcu- 
lations (Bertoglio, Bataille, &: Marion 2001; see also Zank & Matthaeus 
1993), the energy in compressible modes in hydrodynamic turbulence 
scales as ~ M| if Mg < 1. CL03 conjectured that this relation can be ex- 
tended to MHD turbulence if, instead of M^, we use ~ iSV)\/{a^+V^). 
(Hereinafter, we define Va = Bo/y^A-Kp, where Bq is the mean magnetic 
field strength.) However, since the Alfven modes are anisotropic, this 
formula may require an additional factor. The compressible modes arc 
generated inside the so-called Goldreich-Sridhar cone, which takes up 
~ i^y)A/yA of the wave vector space. The ratio of compressible to 
Alfvenic energy inside this cone is the ratio given above. If the generated 
fast modes become isotropic (see below), the diffusion or, "isotropiza- 
tion" of the fast wave energy in the wave vector space increase their 
energy by a factor of ^ Va/{SV)a- This results in 



(SV) 



2 



ml 



Vj + a' {SV)a 



(1) 



where {5V)1^^ and {5V)\ are energy of compressible and Alfven modes, 
respectively. Eq. (1) suggests that the drain of energy from Alfvenic 
modes is marginal along the cascade^ when the amplitudes of pertur- 
bations are weak. Results of calculations shown in Fig. 2 support the 
theoretical predictions. 

We may summarize this issue in the following way. For the incom- 
pressible motions to decay fast, there is no requirement of coupling with 
compressible motions^. The marginal coupling of the compressible and 
incompressible modes allows us to study these modes separately. 



® This reasoning shows that at the energy injection scale when 5B ~ Bo the 

coupling between the modes is appreciable. 

The marginal generation of compressible modes is in agreement with earlier 
studies by Boldyrev, Nordlund, & Padoan (2002) and Porter, Pouquet, & Wood- 
ward (2002), where the velocity was decomposed into a potential component and a 
solenoidal component. A recent study by Vestuto, Ostriker & Stone (2003) is also 
consistent with this conclusion. 

^ The reported (see Mae Low et al. 1998) decay of the total energy of turbulent 
motions Etot follows t^^ which can be understood if we account for the fact that the 
energy is being injected at the scale smaller than the scale of the system. Therefore 
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Figure 2. Mode coupling studies. {a,)left: Square of the r.m.s. velocity of the com- 
pressible modes. We use 144^ grid points. Only Alfven modes are allowed as the 
initial condition. "Pluses" are for low /3 cases (0.02 < (5 < 0.4). "Diamonds" are for 
high 13 cases (1 < /9 < 20). {h)middle: Generation of fast modes. Snapshot is taken 
at t=0.06 from a simulation (with 144'^ grid points) that started off with Alfven 
modes only. Initially, (3 (ratio of gas to magnetic pressure, Pg/Pmag) = 0.2 and Mg 
(sonic Mach number) ~ 1.6. {c)right: Comparison of decay rates. Decay of Alfven 
modes is not much affected by other (slow and fast) modes. We use 216^ grid points. 
Initially, /3 = 0.02 and Ms ~ 4.5 for the solid line and Ads 7 for the dotted line. 
Note that initial data are, in some sense, identical for the solid and the dotted lines. 
The sonic Mach number for the solid line is smaller because we removed fast and 
slow modes from the initial data before the decay simulation. For the dotted line, 
we did not remove any modes from the initial data. Prom CL03. 



3. What are the scalings for velocity and magnetic field? 

Some hints about effects of compressibility can be inferred from the 
GS95 seminal paper. More discussion was presented in Lithwick & Gol- 
dreich (2001), which primary deals with electron density fluctuations in 
the regime of high (3, i.e. (/? = Pgas/Pmag ^ !)• As the incompressible 
regime corresponds to /3 — > oo, so it is natural to expect that for 
/3 » 1 the GS95 picture would persist. Lithwick & Goldreich (2001) 
also speculated that for low [3 plasmas the GS95 scaling of slow modes 
may be applicable. A detailed study of compressible mode scalings is 
given in CL02 and CL03. 

Our considerations above about the mode coupling can guide us in 
the discussion below. Indeed, if Alfven cascade evolves on its own, it is 
natural to assume that slow modes exhibit the GS95 scaling. Indeed, 
slow modes in gas pressure dominated environment (high j3 plasmas) 
arc similar to the pseudo-Alfven modes in incompressible regime (see 
GS95; Lithwick & Goldreich 2001). The latter modes do follow the 
GS95 scaling. In magnetic pressure dominated environments or low (5 
plasmas, slow modes are density perturbations propagating with the 

some energy originally diffuses to larger scales through the inverse cascade. Our 
calculations (Cho & Lazarian, unpublished), stimulated by illuminating discussions 

with Chris McKee, show that if this energy transfer is artificially prevented by 

injecting the energy on the scale of the computational box, the scaling of Etot 

becomes closer to t~^. 
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(b) (d) (f) 

Figure 3. M, - 2.2, Ma ~ 0.7, p ~ 0.2, and 216^ grid points, (a) Spectra of Alfven 
modes follow a Kolmogorov-like power law. (b) Eddy shapes (contours of same 
second-order structure function, SF2) for velocity of Alfven modes shows anisotropy 
similar to the GS95 (ry oc r^/'^ or oc k^/^). The structure functions are measured 
in directions perpendicular or parallel to the local mean magnetic field in real space. 
We obtain real-space velocity and magnetic fields by inverse Fourier transform of the 
projected fields, (c) Spectra of slow modes also follow a Kolmogorov-like power law. 
(d) Slow mode velocity shows anisotropy similar to the GS95. We obtain contours 
of equal SF2 directly in real space without going through the projection method, 
assuming slow mode velocity is nearly parallel to local mean magnetic field in low 
(3 plasmas, (e) Spectra of fast modes are compatible with the IK spectrum, (f ) The 
magnetic of fast modes shows isotropy. Prom CL02 



sound speed a parallel to the mean magnetic field. Those perturba- 
tions are essentially static for a Va- Therefore Alfvenic turbulence 
is expected to mix density perturbations as if they were passive scalar. 
This also induces the GS95 spectrum. 

The fast waves in low (3 regime propagate at Va irrespectively of the 
magnetic field direction. In high (3 regime, the properties of fast modes 
are similar, but the propagation speed is the sound speed a. Thus the 
mixing motions induced by Alfven waves should marginally affect the 
fast wave cascade. It is expected to be analogous to the acoustic wave 
cascade and hence be isotropic. 

Results of numerical calculations from Cho Sz Lazarian (CL03) for 
magnetically dominated media similar to that in molecular clouds are 
shown in Fig. 3. They support theoretical considerations above. 
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Figure 4- left panel: Mach number is 0.35, central panel: Mach number is 2.3 right 
panel: Mach number is 7. The figures are from CL03 



4. What is the scaling of density? 



Density at low Mach numbers follow the GS95 scaling when the driv- 
ing is incompressible (CL03). However, CL03 showed that this scaling 
substantially changes for high Mach numbers. Fig. 4 shows that at 
high Mach numbers density fluctuations get isotropic. Moreover, our 
present studies confirm the CL03 finding that the spectrum of density 
gets substantially flatter than the GS95 one (see also Cho & Lazarian 
2004). Note, that a model of random shocks would produce a spectrum 
steeper than the GS95 one. A possible origin of the flat spectrum is 
the superAlfvenic perturbations created by fast modes within density 
perturbations originated from slow modes. This particular regime is 
clearly identified in a review by Cho, Lazarian & Vishniac (2003) (see 
Fig. 9). It may also be related to the regime of superAlfvenic turbu- 
lence discussed e.g. in Norlund &; Podoan (2003). However, alternative 
explanations of the shallow density fluctuations exist and our current 
work should clarify which process is actually responsible for the unusual 
density scaling that we observe. 



5. How do neutrals affect MHD turbulence? 

An obvious effect of neutrals is that they do not follow magnetic field 
lines and thus produce viscosity. 

In hydrodynamic turbulence viscosity sets a cutoff, with an exponen- 
tial suppression of motion on smaller scales. Below the viscous cutoff 
the kinetic energy contained in a wavenumber band is dissipated at that 
scale, instead of being transferred to smaller scales. This means the end 
of the hydrodynamic cascade, but in MHD turbulence this is not the 
end of magnetic structure evolution. If viscosity is much larger than 
resistivity, there is a broad range of scales where viscosity is important 
but resistivity is not. Over these scales magnetic field will be stretched 
by turbulent motions at larger scales, with the shear from the smallest 
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Figure 5. Viscous damped regime (viscosity > magnetic diffusivity) . Due to large 
viscosity, velocity damps after k ~ 10. (a) Left: Incompressible case with 384^ grid 
points. Magnetic spectra show a shallower slope (Eh{k) oc fc~^) below the veloc- 
ity damping scale. We achieve a very small magnetic diffusivity through the use 
hyper-diffusion. Prom Cho, Lazarian, & Vishniac (2002b). (b) Right: Compressible 
case with 216^ grid points. Magnetic and density spectra show structures below 
the velocity damping scale at fc ~ 10. The structures are less obvious than the 
incompressible case because it is relatively hard to achieve very small magnetic 
diffusivity in the compressible run. From CL03. 



undamped eddies being most important. Indeed, this new regime of 
turbulence has been discovered (see Fig. 5)! 

A theoretical model for this new regime is given in Lazarian, Vish- 
niac, k. Cho (2004; hereafter LVC04). It predicts the spectrum E{k) ~ 
for magnetic field and the spectrum ~ k~'^ for kinetic energy. The 
actual measurements that got somewhat steeper spectra were explained 
in LVC04 as the consequences of a small inertial range available. 

An important prediction in LVC04 is that the intcrmittcncy of mag- 
netic structures increases with the decrease of the scale. This prediction 
was confirmed by numerical simulations in Cho, Lazarian & Vishniac 
(2003b), which showed that the filling factor of magnetic field was 
decreasing with the increase of the wavenumber. 

The effect of neutrals does not amounts only to emergence of the 
viscosity-damped regime of MHD turbulence. Below we describe some 
other effects. 

First of all, it is clear that whether ions and neutrals act as one fluid 
depends on whether the eddy turnover rate t^^^y ~ kvk is longer or 

shorter than the rate t^l of neutral-ion collisions. If t~^y > t~l , neu- 
trals decouple from ions and develop hydrodynamic Kolmogorov-type 
cascade. Indeed, the damping rate for those hydrodynamic motions 
t~l and below the decoupling scale the hydrodynamic motions evolve 
without much hindrance from magnetic field. Magnetic fields with the 
entrained ions develop the viscosity-damped MHD cascade until ion- 
neutral collisional rate gets longer than the dynamical rate of the 
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intermittent magnetic structures. After that the turbulence reverts to 
its normal MHD cascade which involves only ions. 

If t~J^y < up to the scale at which neutral viscosity damps 
turbulent motions, the viscosity-damped regime emerges at the scale 
where kinetic energy associated with turbulent eddies is dissipated. 
Similarly to the earlier case when the when ion-neutral collisions get 
insufficient to preserve pressure confinement of the small scale magnetic 
filaments, outbursts of ordinary ionic MHD turbulence will take place. 
The turbulence will be intermittent both in time and space because of 
the disparity of time scales at which turbTilcncc evolves in the viscosity- 
damped and free ionic MHD regimes. Those predictions should be 
tested with a two fiuid MHD code. 



6. What is the intermittency of MHD turbulence? 

Power spectra do not uniquely characterize turbulence. Very different 
random processes may have the same power spectra. One way to break 
this degeneracy is to use higher order structure functions. 

The p-th order (longitudinal) velocity structure function SFp and 
scaling exponents (^(p) are defined as 

SFp{r) = (I [v(x + r) - v(x)] • r\P) cx r^^P\ (2) 

where the angle brackets denote averaging over x. 

The scaling relations suggested by She & Leveque (1994) related 
({p) to the scaling of the velocity vi ~ l^^^ , the energy cascade rate 
tf^ ~ and the co-dimension of the dissipative structures C: 

Cip) = ^(1 - a;) + C (l - (1 - x/Cr/s) . (3) 

For incompressible turbulence these parameters are g = 3, x = 2/S, 
and C = 2, implying that dissipation happens over ID structures (e.g. 
vortices) . 

Muller Sz Biskamp (2000) applied this model to MHD turbulence and 
attracted the attention of the MHD researchers to this tool. CLV03 dis- 
covered that the magnetic field and velocity has different C: (^'iiagnctic ^ 
^velocity ^ which means that magnetic field is more intermittent than 
velocity. In addition, in CLV03 we found that the value of C changes 
depending whether we use local frame of reference which is perpendicu- 
lar to magnetic field lines or a global one related to the large scale mean 
magnetic field. These complications make us wonder whether the pa- 
rameters of the Shc-Leveque model have the same clear interpretation 
in MHD case as they have in hydro. 
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Figure 6. left pane/:Normalized structure function exponents in perpendicular di- 
rections in the local frame. The velocity exponents show a scaling similar to 
the She-Lcvoquo model. The magnetic field shows a different scaling, central 
pane/: Normalized structure function exponents in the global frame. Note that the 
result for is very similar to the Miiller-Biskamp model, right panel: Magnetic 
structure function exponents, C(p)) in the local frame (not normalized). The observed 
scaling exponents are at least close to the expected asymptote = 0. Figures are 
from CLV03 



7. Summary 

1. MHD turbulence is not a mess. Scaling relations for its modes 
have been established recently. 

2. Fast decay of MHD turbulence is not due to strong coupling of 
compressible and incompressible motions. The transfer of energy from 
Alfven to compressible modes is small. The Alfven mode develops on 
its own and decays fast. 

3. MHD turbulence does not vanish at the viscous scale, provided 
that the fluid viscosity is much larger than resistivity. Instead, a new 
regime of turbulence in which magnetic energy cascade is driven by 
larger scale turbulent motions develops. 

4. Density fluctuations follow the scaling of Alfvenic part of the cas- 
cade only at small Mach numbers. At large Mach numbers the density 
field gets isotropic and has a shallow spectrum. 

5. The intermittency of magnetic field is larger than that of the 
velocity field. The physical meaning of the She-Leveque dimensions is 
not straightforward as it does depend on whether the measurements 
are done in local or global reference frame. 

Acknowledgments We acknowledge NSF grant AST 0307869 and 
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